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MOMENTS OF RANDOM MULTIPLICATIVE FUNCTIONS AND 
TRUNCATED CHARACTERISTIC POLYNOMIALS 

WINSTON HEAP AND SOFIA LINDQVIST 


Abstract. We give an asymptotic formula for the 2kth moment of a sum of 
multiplicative Steinhaus variables. This was recently computed independently by 
Harper, Nikeghbali and Radziwill. We also compute the 2fcth moment of a trun¬ 
cated characteristic polynomial of a unitary matrix. This provides an asymptotic 
equivalence with the moments of Steinhaus variables. Similar results for multi¬ 
plicative Rademacher variables are given. 


1. Introduction 

In the study of the Riemann zeta function there are two probabilistic heuristics 
which have had significant recent attention. One of these is the use of random 
multiplicative functions in problems of an arithmetic nature and the other is the use 
of random matrix theory to model various statistics of the zeta function. 

The study of random multiplicative functions was initiated by Wintner [23] when 
he modelled the Mobius function as the multiplicative extension to the squarefree 
integers of the random variables e p , each of which takes the values {±1} with equal 
probability. This provided a model for the reciprocal of the Riemann zeta function 
and hence an appropriate 1 probabilistic interpretation of the Riemann hypothesis. 
More recently, random models have been used in association with Dirichlet characters 
in the work of Granville and Soundararajan (e.g. in [9, 10]) and also for the quantities 
p lt when p ranges over the set of primes [11, 18, 19]. 

The connection between the Riemann zeta function and random matrix theory is 
well known and has been extensively studied. One of the more remarkable predictions 
of random matrix theory is the Keating-Snaith conjecture [16] regarding the moments 
of the zeta function. This states that 

(1) 1 J |C(| + it)\ 2k dt ~ a(k)g(k)(logTf 
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x We say ‘appropriate’ here since previous models simply used a random ±1 as the coefficients 
which, as objected to by Levy [21], did not take into account the multiplicative nature of the 
problem. 
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where g(k) is a certain geometric factor involving the Barnes G'-function and 


( 2 ) 


a(A:) = n( 1 --) 


P 


40 


rn\ 2 


m =0 


P n 


with dk(n) being the Avfold divisor function. In essence, the reasoning behind the 
Keating-Snaith conjecture can be stated as follows. Since the zeros of the zeta 
function are conjectured to share the same distribution as eigenvalues of a random 
matrix in the CUE, it is reasonable to expect that the characteristic polynomial of 
a matrix provides a good model to the zeta function in the mean. Thus, for an 
appropriate choice of N one could expect 

(3) i f T |<(± + it)\ 2t dt ~ [|A(M, z)\ 2k 

where A (M, z) denotes the characteristic polynomial of the matrix M and the ex¬ 
pectation is over all unitary matrices of size N with respect to the Haar measure. 

Recently, Conrey and Gamburd [7] showed that the asymptotic in (3) holds if 
one both truncates the characteristic polynomial and replaces the zeta function by 
a Dirichlet polynomial of length x = o(T 1 A). This allowed them to deduce that 


lim — 

T->oo T 


E 

n<x 


n 


-1/2 -it 


2k 


dt r^j a(k)c(k) (log x) 


k 2 


k e N 


where a(k) is given by (2) and c(k) is the volume of a particular polytope in M fc h 
This result may be stated in the more general framework of random multiplicative 
functions as follows. 

Given the set of primes, associate a set of i.i.d. random variables {X p }, equidis¬ 
tributed on the unit circle with variance 1. We extend these to the positive integers 
by requiring that X n is multiplicative; that is, if n = rip po(p) then X n = u P xr'. 
We let E[-] denote the expectation. We refer to the X n as multiplicative Stcinhaus 
variables. The association p lt yy X p is then seen to be more than just formal in light 
of the identity 


E 




L n<x 


2q 


= lim — 

I->oo T 



n<x 



which holds for all a G R and q > 0. This can be proved by first demonstrating it 
for q G N and then applying the Weierstrass approximation theorem to the function 

/ : y*-> y q/2 - 

Our main aim is to extend the results of Conrey-Gamburd to more general cr, 
in particular to cr = 0, and to exhibit the connection between moments of random 
multiplicative functions and random matrix theory. 










MOMENTS OF RANDOM MULTIPLICATIVE FUNCTIONS 


3 


Theorem 1. For fixed fceN and 0 < tr < 1/2 we have 


( 4 ) 


E 


n<x 


<k)/3(k) T(2k — 1) fc (1 _ 2o .) , , (fc _ 1)2 

(1 — 2er) 2fc_1 r(/e ) 2 X 1 ° gXj 


where, a(k ) is given fey (2), /3(1) = 1, and 

1 


(5) /3(fc) 


f .. 

. [ 

- k 2k 

TT 17 1 

' i b 2k) 

J(b 2 ) 

11 Si + Sj 
^ i—2 j=k -\-1 ^ 


/c 2k j 

T 2+ - +s ^Y[dsj JJ ^4 

j =2 j=fc+l 


/or k > 2. Here, and bj > 0 /or all j. 

In the significant case a = 0, Theorem 1 has been proved independently by Harper, 
Nikeghbali and Radziwill [13]. It is of interest to note that the constant in their result 
involves the volume of the Birkhoff polytope £> fc . By comparing coefficients we get 
that vol(£>/,-) = k k ~ l fi{k). It is an open problem to determine a closed form for the 
volume of the Birkhoff polytope [22] and a representation in terms of such contour 
integrals may be new. A direct proof of the equation vol(£>/.) = /A _1 /3(fc) can be found 
by applying the methods of section 4.2 to the formula for the Ehrhart polynomial 
given in [4], Also, our methods work equally well in the case <7 = 1/2 and thus by 
comparing coefficients with Theorem 2 of Conrey-Gamburd [7], we get a contour 
integral representation for their constant (see equation (21)). 

It should be noted that the expectation on the line <7 = 0 counts the number of 
solutions (uij ) G N 2k to the equation m\m 2 ■ ■ ■ nik = nrik+i ■ ■ ■ iri 2 k with the restriction 
1 ^ rrij ^ x. In the case k = 2, Ayyad, Cochrane and Zheng [2] computed this 
quantity to a high accuracy. Theorem 1 therefore extends these results to k ^ 3, 
although we do not achieve their level of accuracy. By including the extra condition 
(■ mj,q ) = 1 in the equation, a slight modification of our methods give the following 
asymptotic formula for Dirichlet character sums. 


Theorem 2. Let x fee a primitive Dirichlet character modulo q and suppose q has a 
bounded number of prime factors. Then for fixed fceN, 

( 6 ) 


1 


| 5 >( n ) 

X^XO n<x 



-1 


r (2fc — 1) k 

r(fc ) 2 


(loga:) (fc - 1)2 


as x, q —y oo with x k < q where yo denotes the principal character. 


Our results on the random matrix theory side are as follows. Let U(N) denote 
the group of unitary matrices of size N. For a matrix M G U(N) with eigenvalues 














( e l6j )f =1 let 
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N N 

K{z) = A (M,z) = det(/ - zM ) = JJ(1 - e^z) = ^c M (n)(-z) n 

j= 1 n=0 

where the cm(ji) are the secular coefficients. For N' ^ N we consider the truncated 
characteristic polynomial given by 

A A' (z) = ^ c M (n)(-z) n . 

n<N' 

Let Ej/(jv) [ ■ ] denote the expectation over U(N ) with respect to Haar measure. 

Theorem 3. Let k 6 N be fixed and suppose \z\ > 1. Then for all N ^ k log a; we 
have 

(7) Eu W [|A,„ gI MH ~ (1 _ 

where /3(k) is given by (5) and 

(8) F k (z) =2 F \(1 — k, 1 - t, 2 - 2fc; 1 « |^| -2 ) 
with 2 -Fl 6emy Gauss ’ hypergeometric function. 


One may notice a certain similarity between Theorems 1 and 3. Indeed, by in¬ 
cluding the work of [7] in the case cr = 1/2 we have the following. 


Corollary 1. Let k e N be fixed and let z a be any comlpex number such that 
e i/ 2 -a' Then for 0 < cr < 1/2 and N ^ klogx we have 

2k 


Zrr = 


(9) 


E 


J2 X n/n° 


L n<x 


a{k)c a {k) E(7(JV) [| Aloga: (^cr) | 2fc ] 


where a(k ) is given by (2) and 


C a (k) 


1-e 2 


l-2cr 

1, 


2k— 1 




0 < a < 1/2 
a = 1/2. 


A problem which has garnered some attention recently is to determine the first 
moment of | 'Yh n<x X n \. A conjecture of Helson [14] states that this is o{y/x), but 
this seems doubtful now given the evidence in [5, 13]. Another motivation for the 
present article was to provide a conjecture for the first moment via Corollary 1. 

Let us then assume that Corollary 1 holds for 0 ^ k < 1. Then (presumably) the 
average on the right side of (9) can be taken over matrices of size N = logo; which 
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leads to a computation of the full characteristic polynomial. By an application of 
Szego’s Theorem, Chris Hughes has shown ([15], formula (3.177)) that for \z\ < 1 

do) 

as N —» oo. On applying the functional equation 

A(M, z) = det M(-z) n A(M\ 1 /z) 
with \z\ = e 0 ^ 1 / 2 we obtain the following conjecture. 


Conjecture 1. For 0 ^ k < 1 and 0^cr<l/2u;e have 


( 11 ) 


E 


J2 X n/ n ° 


2k 


n<x 


a(k)F k (e 


1 /2 —cr \ —1 


(1 -e 2CT - 1 )M 2 (l -2 a) 


2k—1 


X 


k(l—2a) 


For k = 1/2 and a = 0 we can compute the constants to a reasonable accuracy. 
The arithmetic factor a(k) admits a continuation to real values of k via the formula 


d k (p m ) 


k + m — 1 
m 


T(k + m ) 
m\T(k) 


We then find that a(l/2) = 0.98849.... The other constants are given by 


Fi /2 (e 1/2 ) 1 = 2 dMi 1; 1 - e x ) 1 = a,gm(l - yjl- 1 + yj 1 - i) 
where agm(x, y) is Gauss’ arithmetic-geometric mean and 

1/4 

= 1.21250... . 


= 0.79099. 


e — 1 


Thus, on combining the constants we acquire the conjecture 


( 12 ) 


E 




L n<x 


0.8769...-v/r 


One can instead consider multiplicative Rademacher variables. In this case, asso¬ 
ciate a set of i.i.d. random variables {Y p }, which are ±1 with uniform probability, 
to the set of primes. Extend this to all positive integers by requiring Y n to be multi¬ 
plicative and non-zero only on the square free integers; that is, Y n = |/r(n)| Y\ p \ n Y p . 

Let 


b(k) = n 


\ k(2k— 1 ) 

p) 



l 

pi 


(13) 
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and 

(14) 


7 {k) = 


07 


2k 


(b 2 k) 


'(Pi) 


n 


i 


2k 


U^ d s, 


l<i<j<2k 




j = 1 


where bj > 0 for all j. We then have the following result. 
Theorem 4. For fixed k e N, k ^ 2 we have 

7 (k)b(k) 2 2k x k (\ogx) 2k *- 3k . 


(15) 


E 


£u 

n<x 


2k 


Let S0(2N) denote the group of orthogonal 2N x2N matrices with determinant 1, 
and let Eso( 2 iv)[-] denote the expectation over S0(2N ) with respect to Haar measure. 

Theorem 5. Let k E N be fixed and suppose z E K, \z\ > 1. Then for all N ^ k logx 
we have 

l(k) 


(16) 


E50(2A) [| Ai oga; (^) | 2fc ] 


(i-lOO 


X 


2k log \z\ 


(log x) 


2k 2 -3k 


where 7(k) is given by (14). 

Corollary 2. For fixed k G N, k ^ 2 and all N ^ k logx we have 


E 




L n<x 


2k 


b{k)2 2k (l - e- l ' 2 Y k ^ so{ 2 N) [|Aio ga: (e 1/2 )H 


where b(k ) is the arithmetic factor given by (13). 

Similarly to the case of Steinhaus variables, we expect that the 1st moment is 
~ Cy/x for some constant c. Unfortunately we have not been able to find an analogue 
of (10) for the special orthogonal group and so cannot make a precise conjecture. 
For some recent results on the order of Y2 n <x^n see [12, 20]. 

2. Asymptotics for Steinhaus variables: Proof of Theorem 1 
2.1. A contour integral representation for the expectation. We have 


E 


J2 x -/ n ° 


2k 


L n<x 


£ 


(rH---n 2k y 


n k+l‘" n 2k 
7lj<X 

We invoke the condition rij < x in each j by using the contour integral 

eds 


(17) 


1 

2ni 


/ y — = 

'(b) s 


1 , y > 1 
0 , y < 1 ’ 


(b> 0 ) 
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with y = x/rij. For each j we take a specific line of integration bj. For reasons that 
will become clear we take bi = e < 1 — 2<j if cr < 1/2 and 5 1 = 2ifcr = l/2. In both 
cases we may take the other lines to be sufficiently large so as to guarantee absolute 
convergence; bj = 2 say (j = 2,, 2k). This gives 


E 


E *•/»'“ = E 


n<x 


1 1 [ .../■ TT f^_\ Sj — 

(n i • • • n 2k y (2m) 2k J [b2k) ] {bl) A1 \ nj J Sj 


n i 

nfc+i---n2fc 


(27 ri] 


2k 


n n 2k 

/ • • • / F k (a + si, ■ • • , cr + s 2k ) TT x 

J{b 2k) J{bi) j=1 


Si ds J 


where 


F k (z 1 ) • •' j k) — 


_ ... 

m-n k = u l u 2k 
n k+l'" n 2k 


Since the condition ri\ ■ ■ -n k = n k +1 • • • n 2k is multiplicative we may express F k {z) as 
an Euler product: 


F k (z 1 ,...,z 2k ) =Yl 


-bm 2fe z 2fc 


p miH-b m k 

=m k+1 +-+m 2k 


P' 


(18) 


k 2k ^ 

j. j. ( ^ J* j^Zi + Zj ^ ^ ^ ^ oj z il+ z j 


pZi + Zj 

y i= 1 j=k+ 1^ 

k 2k 

~-A k (zi ,... ,z 2 fc) n n C(^i + 

i=l j=fc+l 


+^*2 "^.72 


where 


(19) A k (zi, ..., ^ 2 fc) - J i ]JII 


2k 


p L i= 1 j=fc+l 



1 

pmiziH-hm 2t z 2t ' 

miH-b m k y 

=m k+ iH-bm 2 fc 



Upon expanding the inner products and sum whilst referring to the middle line 
of (18), we see that A k (zi ,..., z 2k ) is an absolutely convergent product provided 
+ Zj) > 1/2 for 1 < i < k, k + 1 < j < 2k. 
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We now have 


(20) E 


J2 x -/ n ° 


2k 


L n<x 


(2vr i) 


2k 


' (p2k) 


'(M 


A k {cr + si,..., a + S2k) 


2k 


2k 


X 


n n c(2cr+s,+ xt 

i= 1 j=k +1 j =1 






2.2. The case cr = 1/2. Although the case a — 1/2 has already been investigated 
by Conrey-Gamburd [7], we will go over it with a proof that is instructive for the 
case 0 < cr < 1/2. 

Set a = 1/2 in (20). We write the resulting integral as 


{2'Kl) 2k 

where £ = log a; and 


f f fc 2fc 

/ ••• / 5fc(s 1 ,...,s 2 fc) J] II 

am) Am i=i , =fc+ 


2k 

A—n 

=1 j=k+1 + A j=1 


„Csi dsj 


2k 


Bk(s\, • • • , S 2 k) — Afc(| + Si, ■ ■ ■ , \ + S 2 k) m m (Si + Sj')C(l + s i + s j )• 

i =1 j=/c+l 

This function is holomorphic in a neighbourhood of (0,0,..., 0) and the constant 
term in its Taylor expansion about this point is given by A k (\i ..., |). 

We now make the substitution Sj ha Sj/C in each variable to give an integral of 
the form 


C k 


(27 n. 


2k 


p p k 2k 

/ ■■■ B k (s\/C,.. .,s 2 k/£) nn 

•i(c2fc) Am t=i j=k+ 


1 2k 

—n 

=i j= k+i Si + A j= i 


. dsj 


First, note that we may shift the contours so as to be independent of £, to 9?(sj) = 2 
say. We now truncate the integrals at height T = o(£) and take a Taylor approxi¬ 
mation to Bk(s) about the point (0, 0,..., 0). Then upon letting £ —>■ oo we see that 
this integral is asymptotic to 

fc 2 „ /» k 2k ^ 2k 




£ 


(2ni 


;\2k 


( b 2 k) 


/ n n yryll 

i =1 j=k +1 3 j =1 


rf A 


A short calculation gives A k (\ 1 ..., |) = a(k) where a(k) is given by (2). The re¬ 
maining constant is given by 

i „ k 2k i 2k 


( 21 ) 


a(k) := 


(2 m) 


2k 


ip2k) 


I. II II , - II 

l” 1 ' i=l j=fc+l ^ j=l 
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We may express a{k) as a volume integral and hence recover the constant of 
Theorem 2 in [7]. This is achieved by first writing (s* + Sj) -1 = / 0 °° e~ Xi ^ Si+s ^dxij 
for each term in the product over i,j so that the full product is then given by a 
Ar-fold integral. Upon exchanging the orders of integration and applying (17) the 
result follows. 


2.3. The case 0 < a < 1/2. Returning to our expression for the expectation given 
in (20), we first make the substitutions s 3 H» Sj + 1 — 2cr for k + 1 < j < 2k. This 
gives 


(22) E 


£ 

L n<x 


X n /n° 


2k 


= x k{ 1-2*) 





x A k (a + si,..., a + s k , 1 - a T s k+1 , ■■■ , 1 - a + s 2 k) 


2k 


X 


n n c(i+ Si + sj) 


x° 


dsj 


2k 


i =1 j=k +1 


3 = 1 


3 


n 




ds 3 


j=k +1 


So T 1 — 2 ct 


In the case cr = 1/2, the leading order term was essentially given by the poles at 
Sj = 0. In the present case we must first make the appropriate substitutions to bring 
the leading order contributions to Sj — 0. Only then can we make the substitution 

s j ^ s j/£- 

We first extract the polar behaviour of the integrand. Write 


(23) G k ,a(s i,..., s 2 k) — A k (cr + Si,..., a + s k , 1 — cr + s^+i,..., 1 — cr + s 2k ) 

k 2k 

X II II(* + S i)C( 1 + Si + Sj) 

i—1 j=k+ 1 

so that our integral becomes 


x 


k{l—2cr) 


(2tt7 


2k 


r. r. k 2k 

/ Gk,a(si, ■ ■ ■, s 2k ) ]~[ 

hc*k) J{bi) i= i ,- =fc+ 


~\~ Sj 


■X 


2—1 j=k -\-1 L ' ^ 

k 


x e 


£(siH-l~ s 2fc) 


n 


2k 


dSj 


. , s,- s, + 1 — 2cr 

3=1 J j=k -\-1 J 


The function Gk,a(si ,..., S 2 fc) is analytic in the region 3?(sf + Sj) > —1/2 for 1 < i < 
k, k + 1 < j < 2k. 

We now make the substitutions s 3 n- s 3 — si for k + 1 < j < 2k and s* H> s, + si 
for 2 < z < k. This gives an integral of the form 
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X 


k(l-2a) 


X 


(27ri) 2fc 
r k 2k 

n n 

i=2 j=k -\-1 


\d 2 k) J (&i) 


Gk,a( s I; ^2 + Si ■ ■ • , Sfc + Si, Sfc + i — Si, ... , S2k ~ Si) X 


Si + S j 


, k 2k 

0 C(S2~\ \~S2 k) nA?r n 


els. 


Sl j=2 s i + Sl jtk +1 Sj ( Sj ~ Sl + 1 “ 2cr ) 


Now, for j = 2, 3,... ,2k we let Sj ha Sj/£. This gives the integral 


x k{l-2rr) £(fe-l) 2 


(2vrf; 


2fc 


( e 2k) 


'(bi) 


Gk,a( s 1 , Sl + S2/Si + Sfc/£, — Si + Sfc+i/ — Sl + S2fc/ £) X 


r k 


X 


2fc 


nn 


Sj 

i=2 j=k +1 1 3 


d S2~\ - \-$2k 


ds\ 

Sl 


ds 


2k 


nn?- n 


dsi 


7 2 £ ^ i=fe+i ^ ^ £ 


s j(~r — Si + 1 — 2cr) 


Once again, we may shift the lines of integration in the integrals over S2, S3,..., S2k 
so as to be independent of £; back to 9?(sj) = 2 say, and truncate the integrals at 
some height T = o(£). From the definition of Gk, a given in (23), we see that 

hill Gk,a(s 1, Si + S2I £,..., Si + Sk /£, — Si + Sk+ l/£, • • • , — Si + S2k /£) 

> 00 

=A k (a + si,..., a + si, 1 - a - si,..., 1 - a - si) 

=A fc (0,...,0,l,...,l) 

=Ak(b •••>§) 

where in the last two lines we have used the symmetry of A^. As previously claimed, 
this last quantity is given by (2). The other limits are easily evaluated. 

Thus, as £ A 00 we have 


(24) E 


J2 X n/ U ° 


L n<x 


2k 


a(k)x Hl ~ 2a) C (k ~ 1)2 


(27T i) 2k 


k 2k 

nn 

L i=2 j=k -\-1 


“l - Sj 


S2H- \~ s 2k 


' (p2k) ** 0>l) 

ds 1 


Si(l -2 a - si) 


A: 2k , 

II '- ii 


j=2 j=k +1 
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For the integral over si we push the line of integration to the far left encountering 
a pole at Si = 0. The integral over the new line vanishes in the limit and so 

1 f dsi _ 1 d k ~ l A, _ , k 

2m J {bl) sf(l - 2a- si) fe (k - 1)! ds^ 1 V ° Sl 

1 r(fc + (fc-l))/r(fc) 

(k — 1)! (1 — 2cr) 2fc ~ 1 

_T(2k — 1) 1 

T(k) 2 (1-2 a) 2k ~ v 

The remaining integrals are given by {3(k) of equation (5). 

Note that although it appears as if one should be able to make the substitution 
Sj i —y Sj/C directly in (22) without first shifting the variables by Si, this is not the 
case. Upon truncating the integrals at height T = o(£), the largest error terms arise 
from the ((-factors in G kt „ when they are evaluated close to t = 0. For this to occur 
in all terms of the form £(1 + (s* + Sj)/C) and £(1 + Sj/C) for i = 2,..., k, j = k + 
1,..., 2k, one must have £ « —tj ~ 0 for i — 2,..., 2k and j — k + 1,..., 2k. When 
looking at the error arising from cutting some s* at height T, this large contribution 
clearly is excluded, as one has |£| > T 0 for this i. On the other hand, if one makes 
the substitution Sj H> Sj/C directly in (22) and attempts to cut all integrals at height 
T = o(£), a large error arises from £ « —tj for i — 1,..., k and j — k + 1,, 2k. 



3. Character sums: sketch proof of Theorem 2 


We shall only sketch the proof of Theorem 2 since it is very similar to the proof 
of Theorem 1. Recall the orthogonality property of Dirichlet characters: for m,n 
coprime to q 


1 


X(m)x(n) 

x 


1 if m = n mod q, 

0 otherwise. 


This implies that for all x k < q 


1 

<p(q) 


Y I Y x 


2k 


X m<x 


l 

<p(q) 


EE x( m i ■ ■ ■ m k )x(m k+1 ■ ■ ■ m 2k ) 

X rrii^x 


E 1- 

mi"-m k =mk+i"-m 2 k 
rrii^x 
{mi,q )=1 
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On applying the line integral (17) we acquire 

1 


£ 


1 = 


mi •••ra/ e =771/0+1 •••m 2 fc 
rrii^x 
(m.i,q )=1 


(2ni 


02 k 


„ „ 2k 

/ / Hk,q( s l, ■ ■ ■ i S2k) TT 

•/(ftafc) -Am j=i 


X 


Si ds J 


where 


mi—m fc =m fc+ i—m 2 fc 1 2fc 

(mj,q )=1 


Expressing this as an Euler product gives 


Hk,q{s 1, • • • j <S 2 fc) ^ ^ 


p|g v mi-i -b m k = 

m k +iH-b m 2k 


_ P 


,w»isiS-t m 2k S2k 


- \-m 2k S2k 


n e 


p\q v miH-Fm(.= 

m k+1 -\ -b m 2k 


_ P 


misiH- \-m2kS2k 


n £ „ 

p miH-1 -m k = ** 

m k+1 A -b m 2k 

-1 

^4fc(si ■ ■ • , S2k ) C( s i + s j ) 


~ C kj q(s i, . . . , S 2 fc) C('®* £ ®j")j 




say. Here, the function A k (si,... s 2k ) is that of equation (19). Since the number 
of prime factors of q remains fixed, C k , q (si ,..., s 2 *.) is holomorphic in the same 
regions as A k (s i,... ,s 2k )- The arguments of the previous section now follow, with 
the arithmetic constant being given by 


p\g 


£ 


P 


,mE- \-m k 


=a 


will £ 

p\q 


dk(p 


miH-b m k = 

m k+ i-\ - \-ni2 k 

,m\2\ ~ 1 


-1 


n =0 


jr 


1 

v(q) 


J2xo (n) 

n<x 


l 

<p(q) 



n<x 


(n,q )=1 


1 (p{q) 


<p{q)' q 


x + 0(2^)) 


2k 


Now 
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where oj(q) represents the number of distinct prime factors of q. Since we’re assuming 
oj(q) is bounded this last error term is 0(1) as q —* oo. Hence, 

2k—1 


<p(q) 


n<x 


n 


2k 


y(g) 

q 


2k-l 2k 

X 


< 


q 


<p(q) 

q 


X 


Since this is of a lower order than the main term when oj(q) is bounded the result 
follows. 


4. Moments of the truncated characteristic polynomial in the 
UNITARY CASE: PROOF OF THEOREM 3 


4.1. A formula for the expectation. We begin by recalling the definitions. Let 
U(N ) denote the group of unitary matrices of size N. For a matrix M G U(N) with 
eigenvalues (e t6j )f =1 let 

N N 

A(z) = det(J — zM) = J^J(1 — e l9j z) = cm(»)(- z) n . 

j =1 71 =0 

The coefficients cm(u ) are known as the secular coefficients. We have cm( 0) = 1, 
Cm( 1) = Tr (M) and Cm{N) = det(M). In general, note that these coefficients 
are symmetric functions of the eigenvalues. For N' < N, consider the truncated 
characteristic polynomial given by 

A N'(z) = ^2 c M (n)(-z) n . 

n<N' 

We will compute the expectation of this object as the following multiple contour 
integral. 


Proposition 3. Let k G N. Then for all z G C and N > kC we have 


Et/(jv) [|An(^)| 2fc ] - 


(«!••• U2k ) 


—C 


(2m) 


2k 


n?,. nit+ifi - N 2 ^) “ %(i-%) 


2k 

n 


duj 


where the integration is around small circles of radii less than min(|,z| , 1). 


Our plan is to expand |A(z)| 2fc , push the expectation through, and then use the 
results of Diaconis-Gamburd [8] regarding the expectation of products of the coeffi¬ 
cients Cm(j)- To state their result we must first detail some notation. 

For an rn X n matrix A denote the row and column sums by r* and Cj respectively 
and define the vectors 


row (A) = (rr,... ,r m ), 
col(A) = (ci,... ,c n ). 









MOMENTS OF RANDOM MULTIPLICATIVE FUNCTIONS 


14 


Given two partitions /! = (/x l5 ..., /i m ) and fi = (/*!,..., /i n ) we let JV^. denote the 
number of matrices A with row (.A) = /i and col(A) = /i. The notation (l ai 2“ 2 • • •) 
is used to represent a partition with ai parts equal to 1, <22 parts equal to 2 etc. For 
example, (5, 3, 3, 2,1) = (l^^ 2 ^ 1 ). 

Theorem ([8]). Let (a J -)j_ 1; (h,)* =1 be sequences of nonnegative integers. Then for 
N ^ max ( X]j=i Xq=i .jbj) we have 


Et/(JV) CM(j) aj CM(j) bj — Npii 


where (i = (l ai 2“ 2 • • •) and fi = (l fel 2 62 • • •). 

On expanding the polynomial and pushing the expectation through we get 


Eu(jv) [iA^n 

= ^ Ef/(;v) [cm(hi) • • -CM(nfc)c]i7(nfc + i) • ■ • cff(n 2 k)]z ni+ ''' +nk z nk+1+ ''' +n2k 


: E( 7 (at) [\c M (C)\ 2k ] \z\ 2kC 

+ |z| (2fc_2)£ ^2 E u{N) [c M {C) k ^ 1 c M (m)c^(C) k ^ 1 c^(n)]z m z n d-. 


n^jC— 1 


On taking IV ^ the condition of the Theorem is satisfied for all terms in the sum, 
and can thus be applied. 


Consider the first term. Note that we may write 





'fc 

0 = 1 



k 

•i= 1 


( m u)i.j=i GS fc(C 




( m u)i, 3 =i eB fe(G 


where 
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Similarly, for the second term we may write 

\ z \( 2 k- 2 )c J- E, w [c M (£)*- 1 c M (m)ci7(£)' ! - 1 cw(n)] 2 'V 


m^C 

ri^C—1 


xE?=i EiLi mijyEii ELi 


^ ^ ' ~Z^i=l Z^j=l "Hj ~Z^i=l Z^J = 1 "**J 

n<£-l iGC fc (£,m,n) 


where 


r 2 fc 

C k (C,m,n ) = < (iTiij) e Zg, : ^ 

^ 7=1 


A: A; 

£ for 0 ^ ^ fc — 1, mj-j = m; 

j=i i=i 

A; k 

= £ for 0 ^ j ^ k — 1 , = n 

2=1 2=1 

This set is empty unless m — n in which case we can write the second term as a sum 
of |z| 2 . Continuing in this fashion we see that 


(25) 

where 


E C/(IV) [|A/;(z)| 2fc ] - ^ |£ 

( m >j)i,3=i eD h £ ) 


2Eti E-=i^i 


/- ft; ft; x 

£> fc (£) = < (rriij) e zg 0 : ^ £; ^ £ l. 

f j=i i=i J 

We now invoke the conditions m ij ^ £ with the formula 


(26) 


1 

2ni 


u 


m—C 


du 


m £ £ 


'H=e «(l-w) 1 o, m> C 

which follows on expanding (1 — n) _1 as a geometric series. This gives 

If f 


Eu(iv) [|A/:(^)| 2fc ] - ^ \z 


2 Ei=l ELl m ij . 


(2m) 


2k 


X 


l«2fcl= £ 2fc J\u\\=ei 


X u 


mil+mi 2 H- \-m lk -C m 2 i+m 2 2 -\ - \-m 2k -C „,m kl +m k2 -\ -b m kk -C 


1 


U 


U , 


x n 


mn+m.2ll-bm fcl -£ mi2+m22H bm fc2 --C ,,mu-+™2il-1 -m kk -C 


k+1 


U 


k+2 


■ ■ ■ U 


2k 


2k 

n 


dUn 


= 1 U j (1 M j) 


On collecting like powers and computing the geometric series we acquire Proposition 

3. 
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4.2. Asymptotics for the multiple contour integral. Denote the integral in 
Proposition 3 by I. Thus, 

, 1 f f )~ C TT ,hl 

~ J 'J nt, njh +1 (i - u-i M “it 1 - «i) 

where the contours of integration are positively oriented circles of radii \uj\ = £j. 
We choose l/|z| 2 < £i < l/|z|, £j = £ x l \z\~ 2 5 l J C (j = k + 1,..., 2k) and £j = £ x 8 l J C 
(j — 2,... ,k) with 5j < l/l^l^. With these choices of £j the conditions of Proposition 
3 are satisfied, provided \z\ > 1. We will now perform similar manipulations to those 
in section 2.3 

First, let Uj ha u x l \z\~ 2 Uj for j — k + 1,..., 2k. Then 


I 


log \z\ 

(27 n) 2k 


(u{ k u 2 ---u 2k ) 


-c 




n 


duj 


2k 

n 


dUn 


L = 1 Uj(l - Uj) ^ +1 Uj(l - u x \z\~ 2 Uj) 


Now let Uj ha u i Uj for j — 2,... ,k. After a bit of rearranging we have 


/ = 


x 


2k log |z| 


(u 2 ■ ■ -U 2k ) 


-C 


du x 


(2vr i) 


2k 


u 

dui 


n 


2k 

n 

j=k +1 


duj 


j=2 Uj{l - U l U j) Uj( 1 - Uj){ 1 - Uj \z\~ 2 Uj 


Consider the integral with respect to u 2 ,, u 2k : 


J := 


(u 2 ■ ■ ■ u 2k ) 


-c 


(27 Ti] 


2k—1 


nCn“*+i(i - «i%) 4 “/f 1 - “i u i) 


n 


dUn 


3 

2k 

n 

j=k +1 


duj 


u A l -%)(! 


1 1 £ 

The contours of integration are given by circles of radii £j = o ■' , j = 2,..., 2k. We 
may now choose the to be small but independent of £ since this does not alter 
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the value of the integral. Then, writing J in its parametrised form gives 


J = 


(S 2 ■■ ■ 


(2r 


2k—1 


nf= 2 U^ik+ii 1 - (^j) 1/c e i6i+i9 o) (1 - ui6j /C e i0 i) 


n 


dd; 


2 k 

n 

j=k +1 




(1 — — n 1 1 |^| _2 5,^ £ e* 6 'j) 


Performing the substitutions 9j ha 6j/C and writing Cj = log 8j we have 


rnC 




J = 


g — (c2 + i02-\ - \-C2k +i02 k ) 


£(fc-l) 2 

(2vr) 2k - 1 J_ nC ''' Lc nt 2 n?= fc+ i -C(l - e (c i +ie i +c i +^)/£) 

2k do ,■ 


k 

n 

J=2 


d9i 


n 


Ll (£ _ u ^ e (cj+idj)/C^ 11^ £(£ _ e (cj+i0j)/C^ ^ 1 |,?:| _2 e( c :/+*k))/£) 


We now divide by £( A -' -1 ) 2 and take the limit as £ —» oo. By an argument involving 
dominated convergence we may pass the limit through the integral. Then, since 
£(1 — e~ z t c ) ~ z, we acquire 


J 


£(fc-l) 2 


(1 — Wi) fc-1 (l — u 1 1 \z\~ 2 ) k (27r) 2fc 1 J -00 

p~(. c 2+102 4-fC2fc+*02fc) 


2fc 


d6j 


n* I' pw' 


ni= 2 n i= fc+i (-^ - idi - cj - iOj) j=2 j=k+1 

Finally, we let 6j ha —6 r Upon noting that — Cj is postive, we see that this last 
integral is the parametrised form of /3(k) of equation (5). Therefore, 


P(k)x 2fclog|2| £( fc - 1 ) 2 


du\ 


J (1 — m 1 )^(1 — u 1 qz| 2 ) k u\ 


ft remains to show that 

h : 


(27) 


1 

f 1 

dui 1 

f “L 1 

2ni J 

(1 — Ml) fc (l — U , 1 |z 

_2 ) fc 27ri J 

(1 - Mi) fe (Mi - 

z|' 


-2 \k 


dui 


F(2k — 1) 


F,(z) 


(i — |^|- 2 ) 2fc -i r(k) 2 

where 

(28) F k {z) = 2 Fi(l - k, 1 - k-2 - 2k] 1 - M- 2 ) 

and 2 Pi is the usual hypergeometric function. 
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Note that in I\ we are still integrating over a circle of radius e < l/|z| since we 
made no substitutions in the variable u \. Thus, the only contribution is from the 
pole at U\ = \z\~' 2 . Therefore, 


1 d k ~ l f u k ~ l 
(k — 1)! du k ~ x \ (1 — u) k 

1 /h _ 1 \ rl m Jfc-l-m 

j_bA-iy m+mj/rw rw 2m 

r(fc)“\ m ) (i - \z\~ 2 ) k+m r(m + iy 1 

_1_i_ vi-i ) m ( k ~ 1 \ T ( k + m ) [ 1 \ 

(i ~ kl _2 ) fc r (fc) V ™ / r (m +1) - M 2 / 



As a quick aside, we note that 


fc-i 


(30) 


r(fe) 

' ' m — 


Ecu 


m =0 


m 


1\ T(k + m ) 

/ T(m + 1) 


= 2-Pi(l - k,k\ 1; 1/(1 - |z| 2 )) 


and that this can be written in terms of the Legendre polynomials P n {x) via the 
formula (see [3], section 3.2) 

2 -f\(—A, A + 1; 1; z) — P\( 1 — 2 z). 


Continuing with our manipulations, the last line of equation (29) may be rewritten 
as 


12 k 


(31) 


(\z\ 2 - l) 2fc - 


1 Vi 1) m( k - 1 Y^ k ~ 1 ~ m \ i I~i 2 r 
1 F(k) ^ J { m J r(k-m ) 1 11 J 


12 k 


k -1 


Y(2k — 1) (1 — k) m (1 — k) T 


(\z\ 2 — l) 2fc_1 T{k) 2 ^ ml (2 — 2k) j 


E 


-(i - N 2 ) 


2\m 


where (x) m is rising factorial or Pochammer symbol defined by 


= I 1, m = 0, 

1 x(x + 1) • • • (x + m — 1), m > 1. 

This last sum is the hypergeometric function 2 Pi(l — /c, 1 — k] 2 — 2k ; 1 — |z| 2 ). By 
formula (18) in section 2.9 of [3] we have 

2 Fi(1 — a, 1 - £>; 2 - c;E) = (1 - 2) 6_1 2 .Fi(a + 1 - c, 1 - b ; 2 - c; 2/(2 - 1)). 
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On applying this the result follows. 

5. Moments of Rademacher variables: Proof of Theorem 4 

5.1. A contour integral representation for the norm. Note that E [Y ni ■ ■ ■ Y n2k \ = 
1 if rii - - - 7i2k is a square number and n t is square-free for i — 1 ,..., 2k, and otherwise 
it equals zero. Therefore we have 


E 


E y - 

n<x 


2k 


H n i)\---H n 2k)\- 


Ti\---712k square 


As earlier we invoke the condition rij < x in each j by using (17) with y = x/rij. For 
each j we integrate along the lines bj = 2. This gives 


E 


(32) 


where 


E y - 

n<x 


2k 


E 


2k 


ni--n 2k square 

|mOu)I=i 


(27T i 


;\2k 


I -/ IK-)" — 

'(b 2k ) \ n jJ S J 


„8j d s j 


Ob 


2k 


r r 


2k 

/ 

F k {si,.. 

■, ^ n 

7( b 2k) J( b l) 


i=i 

E 

Tii"'Ti 2 k square 

M n i)' 

"\p(nzk)\ 

< • 

n 2k 


Since the condition ri\ ■ ■ ■n 2 k being square is multiplicative we may express F k (z) as 
an Euler product: 

1 


Fk(z 1) • • • j %2k) /* ^ 


p miH- \-m 2 k even 


n e 


p 


[ m\z\ +—m 2 kZ 2 k 


pz 


l<i<j<2k 

--B k (z 1 ,...,z 2 k) C<A + 


' p z n ^ ^*4 ^ ^ 


l<i<j<2k 


where 

B k (zi, • • •, z 2 k) — J [ 


n i 


p Ll<i<j<2k 


pz 


E 


1 


m iH- \-Tn 2 k even 

COmOl 


P 


miZ!-\ - \-m 2k z 2k 
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Similarly as in the proof of Theorem 1 we get that B k (z) is an absolutely convergent 
product provided + zf) > 1/2 for 1 < i < j < 2k. 

Now we make the substitution Sj K> Sj + 1/2 for j = 1,..., 2k in the second line 
of (32) to get 


E 


E y - 

n<x 


2k 


= X 


(2lTl) 


2k 


[ [ Gk(si, ■ ■ ■ , S2k) II 

M) Bb') i<^< 


1 


l<i<j<2k 


Si + S 


2k 

n 

3 3 =1 


dsj 


x 


Sn + 


where we dehne 


Gk(z i, • • •, Z2k) — Bk(si + 1/2,..., S2k + 1/2) II c(1 +* + »,■)(*+ *,-). 

l<i<j<2k 

The function Gk(z) is analytic in the region St(zi + Zj) > — 1/2 for 1 < i < j < 2k, 
and Gk(0, ■ ■ ■ , 0) = B k ( 1/2,..., 1/2). Finally we make the substitution Sj ha Sj/C 
for j = 1 ,..., 2k to get 


2k 2 —3k 


(2m) 


2k 


J ■■■ J G k ( Sl /C,...,S 2k /C) n 


l<i<j<2k 


Si 


2k 

n 

3= 1 


dsi 


s j/C + 


Shift the lines of integration to be independent of C, say back to 9ft(sj) = 2 for 
j = l,...,2fc and truncate each line at height T = o(C). Computing the Taylor 
expansions and then taking the limit as £ —» oo gives 


E 



6(fc)2 2 V£ 2fc2 “ 3 S(fc) 


where b(k) = B k ( 1/2,..., 1/2) is the arithmetic factor given in (13) and y(fc) is the 
integral given by (14). 


6 . Moments of the truncated characteristic polynomial in the 
SPECIAL ORTHOGONAL CASE: PROOF OF THEOREM 5 

6.1. A formula for the expectation. As for the U(N ) case we begin by expressing 
the expectation as a multiple contour integral. 

Proposition 4. Let k G N, C > 1. Then for all z el and N > kC we have 

\—C 2k 


E50(2IV) [A^(^) 2fc ] — 


(wi • ■■u 2k y 


n 


dui 


(2m) 2k J J III < i< j<2k( 1 - z2 UiUj)f^ 1 u j (l-u j ) 

where the integration is around small circles of radii less than min(|z| -1 ,1). 
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To prove this we must use an alternative method to before since there is no anal¬ 
ogous result of Diaconis-Gamburd [8] for the special orthogonal group. Instead, we 
use the following result of Conrey-Farmer-Keating-Rubinstein-Snaith. 


Theorem ([6]). Let d/i denote the Haar measure on SO{2N). Then for m ^ 1 we 
have 


I SO(2N) 


A(wi) • • • A (w m )dn = 




E 

(ruu) 

s{i,—1} 

\j =i / 


n 


I <i<j<m 1 - w i eiw j ej 


To begin with, we use the integral 
1 

2 m 


(33) 

to write 


f e Ys^ff _ I 1 ) Y > 0 
(0) * 1 o, Y < 0 


(c>0) 


M?) = f- [ A (ze-)e c -L 

2m J [c) s 


for each of the factors A c(z)- Pushing through the expectation then gives 


Es , 0(2A) [A£(^) 2fc ] 



where, for reasons that will become apparent, we take c\ > C 2 > • • • > C 2 k > 
max{0, log \z\} and q — q_i > 2 log \z\ for i = 2,..., 2k. Using the theorem with 
Wi = ze~ Si for i — 1,..., 2k then gives 


E 


SO(2N) 


[A £ (*n= e 


12 Nk 


e f 1 ; — 1 f 


(27Tz) 


;\2 k 


(C 2 k) J(ci) 


2k 

n 

, *=i 


z Net e ~NeiSi 


X 


X 


n 


2k 


(C—N)sj 


1 — y-V'-v • >.i*; 

1<2<j<2Ai j =1 




We would like to show that any term with €j = 1 for some j G {1,..., 2k} gives 
zero contribution. To this end, let n = min {j G {1,..., 2 k} : e 3 = 1}. When 
integrating over si,..., s n _i we only need to keep track of the highest power of e Sn . 
For the si integral, write 


1 

2 e —'U Tfj a,' 



z 2m e m(-si+ejSj) 


m> 0 
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for j — 2,... , 2k. Using (33) we see that m < £ in each sum, so the highest possible 
contribution of powers of e Sn is e Csn . Further, as 62 = • • • = e n _ 1 , this integral 
contributes a nonpositive power of e S2 ,..., e Sn_1 . 

Continuing in this fashion and integrating s 2 ,..., s n _i, we deduce that one is left 
with 



multiplied by some power of z, where the additional terms in f e ( n ^ 1 )- Cs " + • • •) are 
all lower powers of e Sn . For the innermost integral, expanding the factors in power 
series gives 


2k 


2iri J (r \ n 1 — z 2 e Sn+t i s :i 

J j=n +1 


,(C-2N)s n ( (n—l)Cs. 


r e (n-l)Cs n + _ ^ 


ds r 


1 

2n i 


2k 


i c n) \ 7 — 


n 


^ 2 _ S n dnSi 


-2 m J-Sn-ejSj)m \ ( C-2N)s n { (n-l)Cs. 


dSr 


( e (n-l)Cs n + _ 


z “e ^ x ~""e v 

\j=n-\-l m> 0 

The highest possible power of e Sn amongst all terms is e ( 2k - n +( n ~ 1 ) c + c -‘ 2N ) s n which 
is negative for N > Ck and £ > 1. By (33) all terms are thus zero. 

This leaves 

, „ „ , 2 k 


EsO( 2A0 [A/;(z) 2fc ] — 


(2vrF 


2k 


n 


1 — z 2 e 


2 c-Si-S, 


n 

3 = 1 


e Csj —. 


' (c 2 k) J (a) i<i<j<2fc 

In order to arrive at the contour integral, we expand each factor of the product as 

1 


1 — z 2 e 


2/p— So— So 


E 


J2rriij e (-Si-Sj)mij 


Separating the integrals and using (33) then gives 
Eso( 2 A) [a c{zf k ] = 




where 


A k (yC) = \ ( rriij ) G Z>\, 


k(2k-l) 


i —1 2k n 

rriij < £, i — 1,.. ■, 2k >. 
3=1 3=i+l ^ 
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As for U(N) we invoke the conditions Y2 m ij — £ with the formula (26), giving 

If f 


E,so(2 A0 [A j c(^) 2?C | — 


]=£ 


A 52i 


<i<j<2k m ij 


(27T i 


1 2k 


X 


mi2+mi3-\ -hmi2k-£ mi 2 +m 2 3H-hm. 2 ,2fc— 

n 1 u 2 ■ ■ ■ u, 


p2fcl= e 2fc •t|m|=ei 
2k 

mi,2t+m 2 ,2kl-Hri2k_i,2k-£ 

2k 


dUj 


n Wj .(i 


On collecting like powers and computing the geometric series we aquire Proposition 
4. 

6.2. Asymptotics for the multiple contour integral. Denote the integral in 
Proposition 4 by I. Again we perform manipulations similar to those in section 4.2. 
Let Uj i —y \z\~ 1 Uj for j = 1,..., 2k. Then 

2k 


I = 


,j.2k log p| 


(U\ ■ ■ ■ U2k) 


-C 


n 


du,j 


(2vr i) 2k J J - u i u j) fJi u A l ~ M lu j) 

1 1C, 

Next, let u,, ha u/ for j = 1,..., 2k. This gives 


I = 


x 2k log p| 


(ui ■ --u^k) 


-1 


(27rf) 2fc J J 
which can be expressed as 

/ 1 


2k 

n 


duj 


ni<,<,-<»(! { u i u jY^ C ) j =1 Cuj{ 1 — \z\ l u 1 J C ) 
f (ui---U2k)~ 1 


j.2k log p|^2fc 2 — 3fc (27rf) 2 ^ 


2k 

n 


duj 


ril<j<j<2fc£(l _ ( u i u jY^ C ) “T U jY - \ Z \ iyl / c ) 


Here, the contours wind around the origin C times. We now choose the radii of the 
contours to be independent of C and then write the integral in parametrised form. 
Upon taking the limit as C —» oo and pushing the limit through the integrals we 
acquire 

~,2 k log \z\ f*2k 2 —3k 

-7 (k) 


(i-M _1 ) 2fc 


where q(fc) is given by (14). 


7. Concluding remarks 

Admittedly, our evidence for Conjecture 1 is rather weak and there is a certain 
level of ambiguity in choosing the size N of the matrices. However, it is interesting 
that for our choice of N = log x the random matrix expectation seems to capture the 
phase change that we expect to see from the expectation of the Steinhaus variables. 
Indeed, if the conjecture holds for k — 1/2, then we can obtain the order of magnitude 
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predicted by the conjecture for 0 ^ k ^ 1 using Hdlders inequality since we know 
the value at k — 1. Also, it seems a little strange, but not impossible, that one could 
obtain more than square-root cancellation in the case k = 1/2 as conjectured by 
Helson. 

Finally, we note the following argument taken from [5] which gives an upper bound 
on the Steinhaus expectation for k — 1/2. Let 0 ^ u,v < 1 and let S x = Yh n < x Xn- 
By the Cauchy-Schwarz inequality we have 


E[\S X \] 2 <E[|(1 - uX 2 )( 1 - vX 3 )S x | 2 ] • E[|(l - uX 2 )( 1 - vX 3 )\~ 2 ] 
[ "' li 1 -E[|(l — uX 2 )(l — vX 3 )S x \ 2 ] 


Now, 


(1 — u 2 )( 1 — v 2 ) 


E[\(1-uX 2 )(1~vX 3 )S x \ 2 ] 


=E 


(35) 


1 - 2 u$t{X 2 ) + u 2 - 2vlSt{X z ) + 4uuK(X 2 )3fJ(X3) - 2 u 2 v^{X 3 ) 


+ V 2 — 2uv 2 $t(x 2 ) + U 2 V 2 ^j l^a; 


I -1 ,2 2,2 22,2 2,22 

■‘X ( 1 — U + U — 3U + g uv — V + v — UV + U V 


In this last line we have expanded the square of S x and used 


E 


[E 

m,n<x 


X am X b n 


V] 1 ~ —X. 
^ ab 


am=bn 

m,n^:X 


For 0 ^ u,v < 1, the minimum of the function 
f(x,y) = 


l _ „,2 _ 2 ,, , 2 „ 1 „, _ 2„,2 


u + u z — + u 2 — -uu 2 + u 2 v 2 


(1 — u 2 )(l — V 2 ) 
is found to be ~ 0.8164965809.... Taking square roots gives 

E[| S x \\ ^ (1 + o(l)) • 0.903...Vx. 

Of course, further optimisations may prove to disprove conjecture (12). 
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